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Q\ [ It is shown that Renyi statistics provides a plausible basis to describe the 

hadron distributions measured in high energy particle interactions. Gen- 
eralized Boltzmann and gamma distributions obtained by maximization of 
Renyi entropy under constraints on Kolmogorov-Nagumo averages are used 
JpH. to describe the hadron transverse momentum and multiplicity spectra corre- 

Q-i! spondingly. 
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A multi-hadron production in high energy hadronic collisions is a com- 
plicated phenomenon generally modeled by a dynamic system of hadronizing 
quarks and gluons with many degrees of freedom and high level of correla- 
tions. Though using QCD one calculates elementary interactions of quarks 
and gluons at microscopic level, a complexity of the system generally doesn't 
allow an extension of these calculations to predict the observed particle 
spectra. Thus, it is tempting to apply statistical methods to understand 
properties of particle production presented by the experimentally measured 
statistical distributions. In this paper we consider two distributions most 
systematically studied in experiments with pp and pp interactions: the inclu- 
sive charged particle transverse momentum spectrum and charged particle 
multiplicity distribution. 

At low collision energy (s) an exponential behavior of the inclusive single 
particle spectra as function of particle's transverse momentum (P t ) has been 
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observed. These spectra are interpreted using a thermodynamic analogy and 
described by a Boltzmann-type distribution [T] 1 



d? 



a 



= A-e- E ^ T , (1) 

y=0 



d 3 p 

where A is a normalization factor, E t = \Jm 2 + is a transverse energy of 
the produced particle with mass m, and T is a characteristic temperature of 
the hadronizing system. What is a genesis of the exponential distribution? 

The least biased method to obtain statistical distributions, which are 
realized in the nature was promoted by E.T.Jaynes as Maximum Entropy 
Principle (MEP) The MEP states that the physical observable has a 
distribution, consistent with given constraints which maximizes the entropy. 
The Boltzmann exponential distribution (0) arises naturally from a maxi- 
mization of Gibbs-Shannon entropy 

S GS = -Y. P nHPn) (2) 

n 

under a constraint on an average energy of produced hadrons. 

The data on high-s interactions show a deviation from the exponential 
towards a power-law behaviour and can be approximated j2] by 
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= A.(l + ^)-\ (3) 

y=0 



At finite values of the parameter k Q represents a power-law distribution, 
while in the limit of k — > oo the expression (J3J) reduces to Boltzmann ex- 
ponent (JTJ). In plasma physics the /t-distribution (jSJ) is frequently used to 
describe an excess of highly energetic particles with respect to that expected 
from an exponential spectrum [3]. The enhanced probability of the high 
energy particle fluctuations (or an appearance of a heavy tail in a distribu- 
tion) could be viewed as a collective phenomenon resulting from a strong 
intrinsic correlation built in the system. The heavy-tailed distribution (JHJ) 
arises in a framework of non-extensive statistical mechanics |lj as a result of 
maximization of Harvda-Charvat-Daroczy-Tsallis (HCDT) entropy jHj 

1 _ y pi 



1 Since the most data on hadron production are available for central rapidities (y = 0) 
we consider the differential cross sections for central rapidity only. 
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under a constraint on an average energy of the particles. Formally, the ex- 
pression (jHJ) with k = (q— is a non-extensive generalization of Boltzmann 
distribution. 

It is interesting to note that the exponential distribution (JIJ transforms 
in to © when a value of parameter T fluctuates according to a gamma 
distribution jB] 

/(r,A,*) = — ^ . (5) 

The gamma distribution represents a sum of fi + 1 independent exponen- 
tially distributed random variables and arises naturally from maximization 
of Gibbs-Shannon entropy under constraints on the mean values of T and 
ln(T). The latter requirement corresponds to a constraint on a geometric 
mean. 

There are a number of similarities between a behaviour of hadron P t spec- 
tra and hadron multiplicity distributions. The multiplicity distributions are 
observed to be significantly broader [over- dispersed) than the Poisson dis- 
tribution expected for uncorrelated systems. The data on charged hadron 
multiplicities are frequently parameterized by the Negative Binomial Distri- 
bution (NBD) [7j. The NBD arises from Poisson distribution with the mean 
value fluctuating according to a gamma distribution. In other words, the 
NBD results from a discrete gamma distribution with each point smeared by 
a correspondent Poisson distribution. Like the heavy tail in hadron P t spec- 
tra, an over- dispersion of hadron multiplicity distribution is believed to be a 
result of strong long-range correlations in the system. It is therefore tempting 
to find a common statistical approach in which both the heavy-tailed and 
over- dispersed distributions arise naturally. Below it is shown that Renyi 
statistics with Kolmogorov-Nagumo averages provides a plausible basis to 
derive such statistical distributions. 



The generalized Kolmogorov-Nagumo (KN) average of a variable x n is 
defined as 

<X>=<p- l C£ P n<PM), (6) 
n 

where is a monotonic function and P n is a probability of n-th state. For 
x n = ln(P n ) the average (jHJ) represents a generalized measure of information 
(entropy). If one requires the generalized entropy to be additive (extensive) 
and < x + c>=< x > +c, where c is a constant, only linear and exponential 
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functions (j){x) satisfy these requirements. It is practical to represent the 
function <p{x) in a generalized form 

p (i-<j)z _ i 

m = — q > ( ? ) 

which reduces to a linear function in the limit q — > 1. The linear function 
leads to the traditional definition of an average and to Gibbs-Shannon en- 
tropy (0), while the generalized function (J7J) leads to a definition of Renyi 
entropy jS] 

For \q — 1| <C 1 the S q can be approximated by non-extensive HCDT en- 
tropy (0J), while for q — > 1 the entropy (jEJ) reduces to Gibbs-Shannon form (J2J). 
The Gibbs-Shannon statistics has no built-in intrinsic correlations, such that 
for two variables x and y holds < x + y >—< x > + < y >, which is not 
necesseraly for the Renyi statistics. Although Renyi entropy provides 

the most general form of information measure, it has found so far a limited 
number of applications. One of the reasons for that is a little understanding 
of physical meaning of the entropic index q. Since for q = 1 the Gibbs- 
Shannon statistics is restored it is reasonable to assume that a deviation of 
the entropic index value from unit is related to intrinsic correlations built 
up in a system. As we will see below, this intuitive conjecture is compatible 
with the hadronic interaction data. 

In order to obtain a generalized Boltzmann distribution arising in the 
Renyi statistics one has to maximize S q under a constraint on a KN-average 
value. This was done in and in application to a hadron P t spectrum can 
be rewritten as 
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a 



d 3 p 



(q-l)\E t 1 

A-{1- /3 + P-e — —y—, q>l (9) 



with a parameter A = (gT) _1 needed to make the energy dimensionless and 
P being Lagrange multiplier. For |(g — l)XE t /P\ 1 the expression © is ap- 
proximated by the power-law distribution © and the Boltzmann exponential 
distribution (^) follows in the limit q — > 1. 

In Fig. ^ the experimental spectra measured at different pp and pp colli- 
sion energies are shown together with the parameterization ©. The experi- 
mental data come from [TT]. In order to make a comparison with measured 
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differential charged particle cross section Ed 3 a/d 3 p = d 2 a/(2irP t dP t dy) at 
y = the right part of JHJ) is multiplied by E t and m in £" t fixed to charged 
pion mass value. For all data sets a good data description in Fig. ^ is ob- 
served. To compare a quality of the data description using the power-law © 
and the generalized Boltzmann @ distributions the residuals of the corre- 
spondent fits to the data are shown in Fig.|2]for two highest available collision 
energies. It is clear that the power-law function has local systematic deficien- 
cies, while the generalized Boltzmann distribution describes the data equally 
well in the whole available range of the transverse momenta. 

In order to derive a generalized gamma distribution one has to maximize 
entropy S q under constraints on KN-mean values of x n and ln(x n ) 

— =< x > , — ^ =< ln(x) > . 

1 — q 1 — q 

(10) 

In our application, x n = N is the number of produced charged hadrons and 
the resulting distribution is read 

P N = A-((\N)-^ q -^ - (3 + [3e^^)-— , q>l, (11) 

For q — > 1 the (|TT|) is reduced to the gamma distribution ©. 

In Fig. El the experimental charged multiplicity spectra measured at dif- 
ferent collision energies ^2] are fit to the parameterization (jll|) . For all 
data sets a good description by the generalized gamma distribution (JTTJ) is 
found. The fits to the data show that the parameter q increases logarith- 
mically with s both for P t and multiplicity spectra. At the same time the 
experiments show a similar logarithmic rise of a particle-particle correlation 
strength with s [13]. This observation provides a support for the conjecture 
that the entropic index q is related to the intrinsic correlations in system. 
Note, that for = the expression (|llj) transforms in to the generalized 
Boltzmann exponential distribution Q. 

In conclusion, using the Maximum Entropy Principle two generalized sta- 
tistical distributions have been derived from a maximization of Renyi en- 
tropy. These distributions have been applied to describe the experimental 
data on particle production in high energy pp and pp interactions. The 
generalized Boltzman distribution turned out to describe the Pt spectra of 
produced charged particles better than the traditionally used power-law dis- 
tribution. The generalized gamma distribution provides a good description 
of the charged particle multiplicity. 
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Figure 1: Measured charged particle P t spectra overlayed with fits to the 
generalized Boltzmann distribution. 
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Figure 2: A relative difference between the measured charged particle Pt 
spectra and the fit results using the generalized Boltzmann (a and b) and 
the power- law HCDT parameterization (c and d). 




Figure 3: Measured charged particle multiplicity spectra overlayed with fits 
to the generalized gamma distribution. 
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